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It is shown that ponderomotive force can be used to purify entangled states. The protocol is 
based on the possibility to exploit such force for a local quantum nondemolition measurement of the 
total excitation number of continuos variable entangled pairs. 
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I. INTRODUCTION 



o 
o 
o 

Entanglement jl[ has been recognized as one of the most puzzling features of Quantum Mechanics, being also 
O ■ the basis of quantum information processing However, to realize protocols like quantum key distribution or 
quantum teleportation || one needs of maximally entangled states which are not easy to generate (due to loss and 
decoherence) . Thus entanglement purification is often required to distill maximum entangled states from several pairs 
£f~) | of partial entangled states by using local quantum operation and classical communications [Q. 

For qubit systems, efficient entanglement purification protocols have been found j|. Recently, quantum information 
protocols have been also extended to continuos variable systems ||. In particular, in Ref. Q an efficient continuos 
variable entanglement purification protocol has been proposed. Here, a physical implementation of this protocol, based 
on ponderomotive effect is studied. The protocol is based on the possibility of a local quantum nondemolition (QND) 
. measurement of the total excitation number of continuos variables. On the other hand the principle underlying 
a method of measuring electromagnetic energy avoiding photon absorption has been known for approximately a 
hundred years: measure the electromagnetic pressure associated with the energy by an experiment like that of Peotr 
N. Lebedev, and from that pressure infer the energy. Hence, we shall consider such ponderomotive meter Q for 
. entanglement purification. 



a: 



II. THE MODEL 



It is well known J§| that two-mode squeezed states can be generated as intracavity stationary state from the broad- 
band squeezed light provided by a Nondegenerate Optical Parametric Amplifier (NOPA) operating below threshold. 
Let us now suppose that the two modes characterize two different cavities. The latter are also placed in two distant 
sides, say A and B, respectively. Then, the steady state of the cavity modes ai (in the A side) and bi (in the B 
side) can be written as S(r)\0) ai \0) bl where S(r) = exp[r(a\b\ — ai&i)] is the squeezing operator with the squeezing 
parameter r ||. Now, let us consider another pair of modes, say a 2 (in the A side) and &2 (in the B side), in the same 
state, i.e. 5(r)|0) a2 |0)b 2 . Furthermore, suppose that the (bi) modes i— 1,2, belong to the same cavity in the A 
(B) side. For instance, they have the same frequency u> a (loi,), but they distinguish for their polarization. Then, the 
state of the whole system is (£> \ip)2, where \ip)i (i = 1, 2) is a two-mode squeezed state 



\i>)i = Vl-A 2 ^ A» Ql ® \n) bi , (1) 

n=Q 

with A = tanhr. The state |V>)i <8> | V^} 2 represents two pairs of partial entangled states. 

Following the avenue sketched in Ref. [|| we choose to make a QND measurement of the total number of photons 
in the A side. To this end we exploit a ponderomotive meter 0. A suspended mirror provides an intuitive realization 
of such meter JTo| , however, a more realistic apparatus can be realized by using for example a cavity made by a 
piezoelectric crystal coated on both sides pH] . For a piezoelectric crystal, the length variation due to radiation 
pressure can be measured by a suitable electric circuit. 

A simple model to describe the coupling of a ponderomotive meter with the A modes, is obtained as straightforward 
generalization of that presented in Refs. |lJ,[L|. Let c the oscillatory mode characterizing the meter and f2 its 
frequency, then the Hamiltonian (in natural unit) governing the measurement dynamics is 



f 



H = UJ a ^2 a i a i + W & E ^ + ^ C 



E a * ai 



(2) 



with g the coupling constant between signal and meter. The interaction part of the Hamiltonian (fourth term in 
Eq.(Q)), simply represents the effect of the radiation pressure force which causes the instantaneous displacement x 
of the meter ]i4j . Here, x = (c + c f )/\/2, p — —i{c — c^)/\/ 7 2 are the meter position and momentum quadratures 
respectively. 

On the other hand, the meter would be connected with a readout apparatus, and this can be modeled as a heath 
bath . In particular, a meter undergoing quantum Brownian motion can be described by the following master 

equationl|] 



p = -i [H, p] - i- [x, {p, p}\ - - [x, [x, p}\ 



(3) 



where 7 is the damping rate and [3 the parameter proportional to the inverse of the temperature. By assuming the 
meter initially in the vacuum, the initial state of the whole system is 

po = \ft>)i(fl>\ <8> |^)2(^| <8> |0) m <0| . (4) 

Since, from Eq.(^), H commutes with the total number operator in the A side, N — Y?n=i a l a ii one can solve Eq.^) 
by introducing the following characteristic function 



(5) 



where pn,m indicate the density matrix elements of radiation fields built up with 'kef and 'bra' having total number of 
photons N and M respectively. Then, Eq. (|^) can be transformed in a partial differential equation for the characteristic 
function (^). The latter, in turn, due to its linearity allows a simple Gaussian solution. Without going down in the 
details of calculations, it is possible to see that, in the limit t ^> 7 -1 , the following steady state is reached 



00 N N 

p- = E E E a - 

N=0 ni=0 mi=0 



2N 



D{— kN) py^ 1 D\— kN) 



|ni> 0l (mi| ® \N -m) a2 (N -mi\ 
(mi\ <S> \N - ni) b2 {N - m x | 

J th 

where D denotes the displacement operator |17|] , and k = g/Cl. Moreover, 

00 



P^ = (1 



(6) 



(7) 



is the thermal state of the meter. Practically, the dynamics leading to the state (||) corresponds to the 'collapse' of 
wavefunction due to the measurement process . Once 7 becomes large (strong meter-environment interaction) , the 
model provides an idealized limit of real measurement process. Notice that by using the quantum Brownian motion 
master equation no Rotating Wave Approximation has been invoked, differently from Ref. jl6| where the quantum 
optical master equation was used in an analog measurement model. 

The radiation fields state corresponding to the meter outcome x is then obtained by projecting (||) onto the meter 
eigenstate \x) m 1 18 1 , that is 



00 N N 

P ir) -^Y. E E (1- A 2 ) 2 A 2W m (x| [d{-kN)p^D\-kN)\ \x) m 

AT=0ni=0mi=0 

\n x ) ai {mx\ <g> |iV-7ii> 03 (JV-mi| <8) |ni) 6l (m x | <g> \N - n 1 ) b2 (AT - m x \ 

where the normalization constant is given by 

00 

AT 1 = J2 ™(A [d{-kN) P ^ ] D\-kN)\ \x) m , 

N=0 



(8) 



(9) 
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and 



(x| \d(-kN)p^D^-kN)] \x) 



e -f> ~ e -f> 



^ 2™r 

v n— 



■ exp 



+ V2kN 



(10) 



with H n (.) the Hermite polynomials. 

The normalization constant (^) depends on the value of the outcome a;, and gives the meter position distribution 
P(x), that is P(x) = M . Let us now assume for the sake of simplicity f3 = oo, i.e. the meter goes in a coherent 
state of amplitude -kJV @. Then, Fig.(@) shows a typical distribution of outcomes. Different peaks correspond to 

different values of N, from the right to the left N = 0, 1 Notice that in order to distinguish all possible values of 

N one should have the minimum amount of displacement (corresponding to k) greater than the width of a coherent 
state, i.e. 1/2. Therefore, the measurement efficiency increases monotonically with the signal-meter correlation. 

III. ENTANGLEMENT QUANTIFICATION 

At the end of the measurement process the partial density operator can be obtained by tracing Eq.(||) over the A 
side 



(r) 

Pb 



(x)=Tr A {p^(x)} 



oo N 



-Af(l-\ 2 ) 2 ]T 53 X 2N m (x\ [D(- K N)p[fDH- K N)\ \x) m IniMnilSliV-m^JV-nil. (11) 

N=Q ni=0 



One can easily see from the above expression (11) that the state |*ni)6i \n2)b 2 is an eigenstate of the reduced density 
operator 



P%\x) (|n 1 ) 6l |n 2 )6 2 ) = K un2 {x) (|ni) i , l |n 2 )i, 2 ) , 



whose eigenvalues are 



KunM =AA(l-A 2 ) 2 A 2 < ni+ " 2 ) m (x| Id^k^p^dH-kN)] \x) 



(12) 



(13) 



Now, the measure of entanglement for a pure bipartite system is just the Von Neumann entropy of either partial 
density operator j2(J. Hence, it can be written as ||2^| 



£(x) 



53 A »l ."2 ( X ) [ l0 S2 A "l ,"2 ( X )] , 



(14) 



^1,^2—0 



where the dependence from x indicates the dependence from the outcome of measurement. 
Furthermore, the quantity 



T(x) 



So 



(15) 



defines the entanglement increase ratio, once £q identifies the initial amount of entanglement. Of course, if T(x) > 1, 
the state after the measurement is more entangled with respect to the initial one. Prior the measurement we have 
A nii „ 2 = (1 - A 2 ) 2 A 2 (™ 1+ " 2 ) and therefore 



So = 2 [cosh 2 r log 2 (cosh 2 r) — sinh 2 rlog 2 (sinh 2 r)] 



(16) 



which is obviously not depending on x. Furthermore, it results equal to the entanglement of a single pair times the 
number of pairs. 

In Fig. (S) the entanglement increasing ratio T is shown as function of the measured value of position quadrature x 
and coupling constant k. We may see that for k = the entanglement remain the same since the meter is practically 
uncorrelated with the signal. On the other hand, the result x = (corresponding to N = 0) never increases the 
entanglement. For k ^ 0, T grows up approximately as x diminishes; however, one should take into account that the 



3 



probability to get large negative values of x decreases, as can be seen from Fig.([l|). Then, it is useful to introduce the 
probability of success of the protocol 



P s = / dxP(x), (17) 
ir>e 

where integration is performed over the set {x G R|T(x) > e}, with e defining the lower limit for the success (e > 1). 
The quantity Ps is plotted in Fig.(||.a) for e = 1 (solid line) and for e = 1.5 (dashed line). One can see that an higher 
increment of the degree of entanglement is achieved by increasing the value of k, but with a lower probability. 
It is also instructive to consider the efficiency of the protocol. It can be defined as 

H=l-i, (18) 



where 

L. dxT(x)P(x) 

(19) 



t _ J r>e dxT(x)P(x) 
Ir>e dx P ( x ) 

represents the average of the entanglement increase ratio weighted on the probability of success. It always results 
T > 1 and 3 < 1. Fig.(^.b) (solid line) shows the efficiency of the protocol. Contrarily to what one can expect for 
a better resolution of the number of photons, the efficiency does not increase monotonically with k. Instead, there 
exist a particular value of k that maximizes the efficiency, and it is slightly greater than the width of the final meter 
state. In fact, for such value of k, a measurement giving a result, say x, different from zero cannot be identified with 
a particular value of N, but it practically corresponds to any N ^ 0, thus improving the efficiency of the purification 
process. As matter of fact we have seen in Fig.(^) that V exceeds 1 for N ^ (essentially x < 0). However, this limit, 
hence the optimal value of k, also depends on the initial amount of entanglement, i.e. on the squeezing parameter r. 

Fig.(||.b) (dashed line) shows the thermal effects on the efficiency of the purification protocol. In such a case the 
distribution P becomes broader, hence, the protocol results less efficient for k maximizing 3 at /3 = oo. However, 
the effect of temperature is somehow equivalent to decreasing the coupling constant k, in fact for k — > it should 
be T(x) — > 1. Hence the thermal effect can be hindered by increasing n. The extension of previous arguments yields 
k slightly greater than as criteria for efficient protocol. Notice that smaller values of (3 could be considered 
obtaining the same behavior 

Finally, the developed theory can be extended to consider more than two radiation field pairs per side. Consequently, 
the efficiency of the protocol can be improved. Furthermore, in this scheme, all the modes in one side can be considered 
belonging to the same cavity, at least in the side where the measurement is performed (A). However, in that case the 
coupling constant becomes mode depending since we have to account for the mode spacing. 



IV. CONCLUSIONS 



In conclusion, the role of a ponderomotive meter in entanglement purification process is studied. It is found that 
there exist optimal values of the signal- meter coupling constant, contrarily from what one can expect by a simple 
analysis on indirect measurement. The purification process results sensible to the thermal effects introduced by the 
meter, however, appropriate signal-meter coupling reduces such detrimental effects. It is worth noting that a wide 
range of n values are accessible for movable mirrors with the present technology p2| . Otherwise, microfabricated 
cantilevers [ p3[ could be also employed. Furthermore, a large number of pairs involved in the process renders it more 
efficient. 

Since radiation pressure allows to couple macroscopic objects (like mirrors) with microscopic systems (like cavity 
modes), it provides a unique bridge between quantum and classical worlds. Hence, optomechanical systems are worth 
studying, also for quantum information processing. 

Finally, ponderomotive effects come out also in the dispersive regime of atom- field interaction p4| ] , where the atomic 
momentum becomes sensible to the photon statistics, and a scheme analogous to that described can be developed 
with atoms crossing a cavity. In such a case the thermal effects should be less significant. 
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FIG. 3. The success probability Ps is shown in figure a) vs the coupling constant k. The solid line is for e = 1 and the 
dashed line for e = 1.5. The efficiency S of the purification protocol is shown in figure b) vs the coupling constant k for e = 1. 
The solid line is for (3 — oo and the dashed line for j3 — 1. In both figures a) and b) r = 0.3. 
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